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A NOTE ON FUNCTIONAL A POSTERIORI ESTIMATES FOR ELLIPTIC 

OPTIMAL CONTROL PROBLEMS 

M. WOLFMAYR 


Abstract. In this work, new theoretical results on functional type a posteriori estimates for 
elliptic Optimal control problems with control constraints are presented. More precisely, we 
derive new, sharp, guaranteed and fully computable lower bounds for the cost functional in 
addition to the already existing upper bounds. Using both, the lower and the upper bounds, we 
arrive at two-sided estimates for the cost functional. We prove that these bounds finally lead to 
sharp, guaranteed and fully computable upper estimates for the discretization error in the state 
and the control of the optimal control problem. First numerical tests are presented confirming 
the efficiency of the a posteriori estimates derived. 


1. Introduction 

During the last couple of decades, the optimization of systems governed by partial differential 
equations (PDEs) has become more and more important in research and application, for which 
Lions has definitely paved the way with his work [53] in 1971. Books considering PDE-constrained 
optimization are, for instance, Hinze et al. na, Trdltzsch (35] , Borzi and Schulz |3| and Leugering 
et al. |22l [5T] . Besides the PDE-constraints, the optimization problems often include control 
constraints given through a non-empty, convex and closed subset of a Hilbert space. In many 
cases, the set of admissible controls is represented in terms of inequality constraints (or box 
constraints) imposed on the controls. 

A posteriori error estimates and adaptive methods for elliptic optimal control problems are the 
topic of many works, see, e.g., niiiMiiiiiniiiniiiiiisniiiij, which are mainly on residual-type a 
posteriori error estimates. Regarding a posteriori error estimates for control constrained optimal 
control problems, we also refer to the recent works [13 HD- There are different approaches to a 
posteriori error estimation. Besides the residual-type estimates, there is also the class of functional 
type a posteriori error estimates. These techniques were introduced by S. Repin in the 90’s, see, 
e.g., [53[311 [31 [31133HD- Later the functional type a posteriori estimates were also considered 
and obtained for optimal control problems, see [710 as well as the books [351 [55] and the references 
therein. The benefit of functional type a posteriori estimates is that - as the name indicates - they 
are only derived by functional methods. Therefore, these estimates do not depend on the mesh 
and provide guaranteed upper bounds for the discretization errors. 

In [7], sharp, guaranteed and fully computable upper bounds (majorants) for cost functionals 
of distributed elliptic optimal control problems were already presented. Now, we want to complete 
the functional type a posteriori error analysis of distributed elliptic optimal control problems by the 
derivation of guaranteed and fully computable lower bounds (minorants) for the cost functionals. 
The optimization problems of this work also include control constraints. The presented minorants 
are not only fully computable but sharp. Moreover, we prove that the new minorants and the 
already obtained majorants of [7] can be used in order to derive functional type a posteriori 
error estimates for the discretization error in the state and the control. The properties of the 
majorants and minorants are transferred to the majorants for the discretization error, i.e., they 
are guaranteed, fully computable and sharp. 

The paper is organized as follows: In Section [3 we discuss the optimal control problem and 
its optimality system as well as present the already known results on majorants for the cost 
functional. Then, Section [3 is devoted to the derivation of new, fully computable minorants 
for the cost functional. Both results, the older ones on majorants as well as the new ones on 
minorants, together lead to guaranteed and fully computable upper bounds for the discretization 
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error in the state and the control. These majorants are presented in Section |4l They attain their 
exact lower bound on the exact solution of the optimal control problem and, hence, are sharp. 
In Section O we present the corresponding results for the optimal control problem in the case 
without any inequality constraints imposed on the control. Section [5] is devoted to the finite 
element discretization of the optimal control problem and its iterative solution method in order 
to derive an approximation of the solution. We present first numerical results in Section [71 and 
finally draw some conclusions in Section |S1 


2. The Elliptic Optimal Control Problem 

Let fl C d G {1,2,3}, be a bounded Lipschitz domain with boundary T := and, let us 
denote the state of our optimal control problem by y and the control by u. Given yd G 
Ud G / G and A G R+, we consider the following distributed elliptic optimal control 

problem: Minimize the cost functional 

(1) J(y(v),v) := ]^\\y - VdW^ + ^\\v - Udf 

over y G V := and v G C/ad C subject to the boundary value problem 

—div(iy(x)Vy(x)) = f(x) + v(x) x G ^2, 

y(x) =0 a; e r. 

The set of admissible controls C/ad is given by 

(3) C/ad = (f S L^(fi) : Ua < V < Ub a.e. in fl}, 

where Ua,Ub G and Ua{x) < Ub(x) for almost all x G ft. The diffusion coefficient p(-) is 

assumed to be measurable, uniformly positive and bounded, i.e., satisfies the assumptions 

(4) 0 < p < v{x) <v, X G fl, 

where p and p are constants. The positive regularization parameter A provides a weighting of the 
cost of the control in the cost functional In these work, we denote by (•,•) and || • || the 

inner products and norms in L‘^{fl), respectively, whereas the standard inner products and norms 
in H^{ft) are denoted by (•,-)i || • ||i, respectively. 

Our goal is to approach to the desired state function yd and to the desired control Ud as close as 
possible by finding a suitable control function u. Note that, in many problems, the desired control 
is given by Ud = 0. The optimal control problem ([I])-© has a unique solution (for the proof, see, 
e.g., [HISE]) that can be also derived via the optimality conditions. Hence, the optimal solution 
of the optimality system is equivalent to the solution of the optimal control problem (|I])-(ini). In 
order to formulate the optimality system, we consider the following Lagrange functional for the 
minimization problem: 

(5) C{y{v),v,p{v)) ■.= J{y{v),v)+ / {d\v{v{x)Vy{x)) + f{x) + v[x))p{x) dx, 

JQ 

where p denotes the Lagrange multiplier (adjoint state). The Lagrange functional ([5|) has to be 
understood in the weak sense. It has a saddle point, see, e.g., IHISHl. Hence, the corresponding 
(optimal) solution {y,u,p) G V x C/ad x V satisfies the system of (first order) necessary optimality 
conditions 


£p{y{u),u,p{u)) = 0, 
Cy{y{u),u,p{u)) = 0, 

{Cu{yiu),u,p{u)),w - u)l2(q) >0 Vui e C/ad, 
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which can be written in the weak form as follows: 


(6) 

<1 

< 

S 

II 

(/ + u)w dx 

Vw 

G 

V, 


Jn 

Jn 





(7) 

<1 

Vwdx = / 

(y - yd)wdx 

V w 

G 

V, 


JQ, 

Jn 





(8) 

/ {p + X{u - Ud)){w 
Jn 

— u)dx >0 


V w 

G 

Uad, 


for (optimal) y,p G V and u € Uad- 

Remark 1. In the unconstrained case, i.e., Uad = the variational inequality 0) of the 

optimality system simplifies to the equation Cu{y{u),u,p{u)) = 0. Hence, the optimality condition 
13) is then simplified to 

(9) p + \{u — Ud) =0 in Vl. 


Majorants for the cost functional. Guaranteed and fully computable upper bounds for the 
cost functional have already been presented in [7], i-e., 

(10) J{y{v),v) < J®{a,/3;r],T,v) VvGUad 


and for arbitrary a, jd > 0, rj G V and 

T G H(div,fl) := {r G [L^(0)]^ : divr G 


Although the admissible set is different in [7] and the elliptic boundary value problem is stated 
without a diffusion parameter u(-), the majorants can be analogously derived for the optimal 
control problem ((IJ-lISl)- Including a diffusion parameter that meets the assumptions (|3]) in the 
model problem ([2]), the majorant J'® is given by 


:=^-^\\r] - ydf 


(1 + a)(l + l3)Cp . 


( 11 ) 


2aiy^ 


(1 + a)(l +/3)C'|. 


div' 


|T-uV?7f 

+ ^h-Ud\\ 


where Cp > 0 is the constant coming from the Friedrichs inequality. The parameters a, /3 > 0 
have been introduced in order to obtain a quadratic functional by applying Young’s inequality. 
The arbitrary functions p G V and v G Uad can be taken as the approximate solutions of the 
optimal control problem ([I|)-(I31) and r G i7(div, fl) represents the image of the exact flux uVr;. 
For the derivation of (EH), the following estimate for the approximation error has been used: 

(12) l|Vy(u) - V 77 II < ^ (||r - uViyH + CpWf + u + divx||). 


The majorant EH) provides a sharp upper bound of the cost functional, if it is minimized over 
p,T,v and a,/3 > 0, i.e., 

(13) inf J®{a,/3;p,T,v) = J{y{u),u), 

rjG H{div, 

V£Uad,Cii,0>O 


since the infimum is attained for the optimal control u, its corresponding state y{u) and its 
exact flux vVy{u), and for a going to zero. Hence, (fT51) states that the exact lower bound of 
the majorant EH coincides with the optimal value of the cost functional of the optimal control 
problem. Therefore, we have the estimate 

(14) J{y{u),u) < J®{a,P;p,T,v) V p G V, t G H{div,D,), v G Uad, a, ft > 0, 


see [TIUTH]. 
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3. Minorants for the Cost Functional 

In this work, we enrich the derivation of guaranteed upper bounds for the discretization er¬ 
ror in the state and the control of problem ((I|)-(1S]) by obtaining fully computable lower bounds 
(minorants) for the cost functional J. For any rj G V, we have that 

J{y{v),v) = ^\\y- v\\^ + yd) dx +^h- 2/df + ^11^ - 

for all V G C/ad- Since — > 0, we can estimate J from below by 

(15) J{y{v),v) > - yd\\^ + ^\\v - Udf + J^{y - v) iv - yd) dx. 

Let Prj G V he the adjoint state corresponding to rj G V. Hence, solves the equation 

(16) / vVpr^-Vwdx= / {r] — yd)wdx MwGV. 

Jn Jn 

By using (1161) . it follows for (fT51) that 

J{y{y),v) > i||?7-2/df + - UdW^ + jy - V y) ■ vVp^dx 

1 X f 

= - ydf + ^\\v - Ud\\^ + / {lyVy- vVrj) -Wpndx. 

J Q 

Since y = y{v) solves the variational formulation 


(17) 


/ vVy-Vwdx= / {f + v)wdx MwGV 

Jq. Jq 


of the boundary value problem Q, we obtain 

J{y{v),v)>^\\r]-ydf + ^\\v - UdW^ + [ {f + v) p^ dx - [ vVri-Vpr,dx. 

For any x G H{div, H), the identity 

(18) / divTwdx = — / T-Vwdx MwGV 

Jn Jn 

is valid, which yields 

J{y{v).v)>^\\'q-yd\\‘^+ ^\\v-Ud\\'^ + f {f + v)pndx- f vVyVpr^dx 
+ / div TPr^dx + / T • S/pr) dx 

Jn Jn 

= 7 ;h-ydW^ + ^\\v-UdW'^ + [ {f+ v + divT)pr,dx+ [ {t - vV-q) -Vprjdx 
^ ^ Jn Jn 

for all V G C/ad- Moreover, 

J{y{u),u)= inf J{y{v),v) >l-\\q - ydW^ + [ {f + div t) p^^ dx + [ {t - i/Vq) ■ Vp^^ dx 

z Jn Jn 


-I- inf 
veiled 


[ vpn dx + -^llw - Md||- 
Jn ^ 


The control v G C/ad can be defined via the projection formula 

IP’[a,h](M) := min{6, max{a. It}} 

for all a,b gM. with a < b and u gM., which projects K on the interval [a, b]. It is given by 
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where p = p{v) is here the adjoint state corresponding to v and y{v). Let us denote by '^Pri S ^ad 
the control corresponding to the adjoint state Prj, i.e., 

( 19 ) Vp^{x) = ~ '^Pr]{x)} 

for almost every x G fl. By adding and subtracting Vp^, we obtain 

1 A A f 

J{y{v),v) > -h-yaf + - “df + - VpJ'^ + X {v - VpJ {vp^ - Ud) dx 

+ / {v-Vp^)ppdx+ / (/ + Up,+divr)p^da;+ / (t - uVij) ■ V Pp dx 
J Cl J Cl J Cl 

for all V G U^d- Since ^ ||?; — Vp^ |p > 0, we get that 

1 A f 

Jiy{v),v) > -Wv-ydf + ^\\vp^ -udf + (v-vp^) {Pp + X{vp^ - Ud)) dx 

+ / {f + vp„ + div t) ppdx+ (t — vVp) ■ Vpp dx. 

J Cl J Cl 

Due to the variational inequality 


( 20 ) 

we obtain 


/ {v-Vp^){pp + X(vp^-Ud)) dx>Q Vue C/ad, 
Jn 


Jiy{v),v)>^\\p-ydf+ ^\\vp^-Udf + [ {f+ Vp^+dwT)ppdx+ [ [t - vVy) -Vppdx. 

^ ^ Jn Jn 

Now, let us introduce an arbitrary ^ e 1/ and its corresponding control V(^ e C/ad, which can be 
computed by the projection formula as follows: 


( 21 ) 


1 


Uf (a;) = - -C(3:)}- 

Hence, the following variational inequality is satisfied: 

(22) / (u - u^) (C + X{v(; - Ud)) dx >0 Vu e C/ad- 
Jn 

Adding and subtracting e C/ad leads to 

1 A /* 

J{y{v),v) >-\\p-ydf + ^\\vc-Udf + >^ {vpr, - ^c) (u^ - Ud) dx 

+ / (/+ + divT)pp d£c + / [t - vWrj)-Vppdx + / {vp^ - v^jppdx 

J Cl J Cl J Cl 

>\h-ydf + “^c) (Pp-Odx 

+ / if + + div T)ppdx + / [t - vS/r]) -Vppdx, 

Jn Jn 

since -IHup^ — > 0 and the variational inequality (1^^ is valid. Next, we add and subtract the 

arbitrary C, (zV leading to the following estimate: 

J{yiv),v) >^h - yd\\^ + ^\\vc - UdW'^ + J^{vp^ -vc) {pp-C)dx 

( 23 ) + / {f + V(^ + divr) ( dx + / {t — v\Jp) ■ C, dx 

Jci J Cl 

+ / {! +VC, + div t) (jpp - C) dx + / {t - vVp) ■V{pp-C)dx. 

J Cl J Cl 

The following result provides an estimate for the error in the control: 
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Lemma 1. Let Vp^ € C/ad o,nd V(^ € C/ad satisfy the variational inequalities and respec¬ 
tively. Then, the error between Vp^ and V(^ can be estimated by 


(24) 


\\Vp„ -l^cll < ^lbr,-CI|- 


Proof. Adding the variational inequalities 

/ -'>^c) iC + - Ud)) dx > 0 and / {vc; - VpJ {ppX{vp^ - Ud)) dx > 0 

Jn JQ 

yields the inequality 

/ (^Pr, - (C -Pv)dx>X (Vp^ - V(;f dx. 

Jfi J fi 

By applying the Cauchy-Schwarz inequality, we obtain 

Ibp, - ^^clllbp - Cll > / K,-^^c) (C-Pp) > A|bp, 

Jn 

which finally leads to the estimate (EH). □ 


In order to formulate, a computable lower bound for the cost functional, we need to prove a 
computable upper bound for the error in the adjoint state, which is presented in the following 
theorem: 


Theorem 1. Let yd € L'^iLl) he given and let pp G V meet equation il6]) with rj G V and v(-) 
satisfying assumption 0- For any C, gV, we have that 

(25) l|V(pr,-C)ll < ^{CpWv-yd + divpW + ||p-j^VC||), 

where p G iC(div, fl) and Cp > Q is the constant coming from the Friedrichs inequality. 


Proof. Since the bilinear form of problem (1161) is elliptic with ellipticity constant ^ and applying 
the Cauchy-Schwarz and Friedrichs inequalities, we get that 


K\\y{pp 


Oil < sup 

0^w£V 


In ~ C) ■ dx fn(d ~ Vdfw — nX/if ■ X/wdx 

llVrcll ^ oSv 


= sup 

O^wGV 


fn(d ~ Vd + div p)w + (p — i^VC) • Vwdx 


/ lb- 2 /d + divp||||u;||-h ||p-:/VC||||Vw|| 

^ sup .. .. 

l|vin|| 

< C'f lb-2/d + div pII + \\p-pX/C\\ 


where p G Lf (div, il) satisfies identity (ITHll . Hence, it follows the estimate 


□ 


By using the Cauchy-Schwarz and Friedrichs inequalities as well as Lemma [TJ we can further 
estimate the inequality (l23l) from below as follows: 

J(y(y),v) > ^lly- Pdf + ^Ihc - + J (f+ v^+ divT)Cdx + J (t -vVrf) - Vfdx 

- Ibp, -^^clllbp - Cll - ll/ + ^^C +divT||||p^ - Cll - I|r - j^Vr7||||V(p^ - Oil 

>-^lb-2/d|P + ^lbc-'“dll^+ / (f + V(^ +div T)Cdx+ [ (t -vV-q) -Vfdx 
^ ^ Jn Jn 

- jWPr, - Cf - (Cf ||/-hwc -f divrll 4- ||r - ^V?7||)|| V(p^ - Oil 

>-^lb-2/d|P + ^lbc-'“dlP+ / (f + V(; + div t) C dx-G [ (t - vVq) -VQdx 
^ Jn Jn 

- ^l|V(p^ - Of - (Cf ||/-4?;c +divT|| -4 ||r - ^V?7||)||V(p^ - Oil- 
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Now, applying Theorem [T] yields the following estimate: 

(26) J{y{v),v) > V?7,C €V,T,p€ H{div,n),v^ = - ^C} 

with the minorant 

C, r, p, = i ||?7 - ^ \\vc - udf + + div r) ( dx 

(27) + [ {t- vS/rii) • VCd® - - {Cp ||r? - i/d + divp|| + \\p- uVC||) 

Jn iL 

X {^h - Vd + di-vp\\ + ^||p-uVC|| +Ci;’||/ + i;^+divr|| + ||r-uV7?||^ , 
where T,p G H{div,il). Note that the minorant is fully computable. 

Theorem 2. The exact upper bound of the minorant defined in i27\) coincides with the optimal 

value of the cost functional of problem or, equivalently, of the optimality system 

i.e., 

(28) sup J^{p,C„T,p,vtf) = J{y{u),u). 

77, G V, T ,pG (div, r 2 ), 

Proof. The estimate (1^5)) is valid for all v € Uad- Hence, also for the exact solution u, i.e., 

(29) J{y{u),u)= inf J{y{v),v)> J^{r],C,,T,p,vcf) 

ad 

for all T],( G V, T, p G i7(div, H) and for the control V(; G Pad, which depends on ( and can be 
computed by the projection formula m- For the exact solution Vi^ = u, y = y(u), C, = p{u), 
T = vVy{u) and p = vVp{u), the estimate is sharp, i.e., 

J^{y{u),p{u),vVy{u),vVp{u),u) = ]^\\y - yaf + ^||u- = J{y{u),u). 

□ 


4. A Posteriori Error Estimates for Control and State 

In this section, we will derive guaranteed upper bounds for the discretization errors of the 
control and the state measured in the following norm: 

(30) \\\u-v\\\^ := i||y(u)-y(u)f+ ^||M-uf, 

making use of the ideas based on the work by Mikhlin [53] but generalized for the class of optimal 
control problems, see also |35] . 

Theorem 3. For any control function v G Pad, we have the estimate 

(31) |||u- u||p < J{y(v),v) - J{y{u),u). 


Proof. We compute the difference 

J{y{v),v) -J{y(u),u) = ^\\y{v) - ydf - ]^\\y{u) - yaf + ^11^- “df “ ^11“ “ «d||^ 


1 

2./n 


[ iy{v) + y{u) - 2yd){y{v) - y{u)) dx + ^ [ {v + u - 2ud)iv - u) dx 
Jn ^ Jn 

= ^ [ iy(.v)-y{u) + 2y{u)-2yd){y{v)-y{u))dx + - [ {v - u + 2u - 2ud){v - u) dx 
^ JQ ^ Jn 

= \\\y{u) - y{v)\\'^ + - yd){y{v) - yiu)) dx + - uf ^ ~ 
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Since the exact adjoint state p{u) S V fulfills (O, we obtain the equation 
J{y{v),v) -J{y{u),u) 

= h\y{u) - yivW + I iy\7p{u){yy{v)-\7y{u))dx +-\\u-vW^ + X I {u - Ud){v - u) dx 
JQ. ^ Jn 

= ^Il2/(w) - + [ Vp{u){vVy{v) - vVy{u)) dx + ^||m - nf + A [ {u - Ud)iv - u) dx. 

JQ ^ Jn 

From equations ([HI) and (fT71) follows that 
J{y{v),v) -J{y{u),u) 

= -^IlyCw) - y(^^)||^ + [ p{u){f + v- f -u)dx + ^\\u- u|p + A [ {u- Ud){v - u) dx 
^ Jn ^ Ja 

= |||m —'cjlP + / {p{u) + X{u — Ud)){v — u) dx. 

Jn 

Since the variational inequality ([5]) is satisfied, i.e., 

/ {p{u) + X{u — Ud)){v — u) dx > 0 VliGUad, 

Jn 


we finally obtain the estimate (ED). 


□ 


Theorem 4. For any C, £ V, let vq € Uad be given by the projection formula i21\} . Then, we 
obtain the following error majorant: 

(32) |||m-'!^cIIP ^ M®{a,l3-,riX,T,p,VQ) := (a,/3; ry, r, {p,C„T, p,V(f) 

with 

M®{a,l3;p,(,T,p,v^) = ^\\v - ydW^ + ~ 


(l + a)(l + /3)C| 
2af}vf 


11/+ + divx|p - / {f + V(; + div t) C dx 

Jn 


- [ (t- vX/y) ■X/C,dx + - [Cf ||?7 - yd + divp|| + ||p- z/VC||) 

Jn E 

X yd + divp|| + ^\\p-vVC,\\ +C'f||/ + ?;c +divT|| + ||t- i/VryH^ 

for arbitrary rj £ V, t, p £ H (div, O) and a,j3>Q. 

Proof. Using ED, we obtain the estimate 

Ill'u-'ycIlP < - J{y{u),u). 

Applying (flUll and (1^91) finally leads to the estimate 


□ 

Proposition 1. The majorant At® (a, /3; ij, f, t, p, V(f) defined in ED attains the exact lower bound 
on the exact solution of the optimal control problem or, equivalently, of the optimality 

system ©-(ID, *.e., 

CPV m M®{a,p-,r],C,T,p,vc)=0. 

77,G V ,T,pG-Ti (div,r 2 ), 

The infimum is attained for V(^ = u, rj = y{u), ( = p(u), t = vVyiu) and p = vVpiu). 

Proof. We have that 


At® (a, /3; y{u),p{u), vVy{u), vVp{u), u) = -\\yiu) - ydf, 


which is zero if we let a go to zero. 


□ 
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Although the majorant M.® is a guaranteed, computable and sharp upper estimate for the 
discretization error in the combined norm, it only decreases with order h when discretizing the 
mesh. However, the combined norm ||| - ||| is an L^-norm, and, hence, decreases with order 
So the majorant M® is an overestimation for the combined norm. Now, we introduce another 
norm which is a weighted 77^-norm for the state depending on the corresponding control. More 
precisely, we derive an estimate for the discretization error measured in the following norm: 

(33) |||u - ulll^ := ^\\y{u) - y{v)\\'^ + ^-||V 2 /(u) - V?/(u)f. 

Theorem 5. For any control function v G Uad, we have the estimate 

o \ 

(34) |||u- ■cllli < Jiy{v),v) -J{y{u),u) + (||r- uV7?|| + CpWf + v + diYrWf . 

Proof. Let (5 > 0 be an arbitrary but fixed parameter. Adding and subtracting V 77 as well as 
applying triangle inequality for ^-j||V?/(m) — V?/(u)||^, we derive the following estimate: 

^||V 2 /(u) - V 2 /(u)f < ^ (||Vy(u) - Viyf + ||Vy(u) - Viyf) • 

Using adding and subtracting v as well as applying twice triangle inequality, we arrive at the 
estimate 

^||V 2 /(u) - V?/(u)f < uViyll +C'F||/ + u + divT|| +C'f||u-z;||)^ 

+ (Ik - + Cull/ + w + divrll)^^ 

- +CF\\f + v + divx||)2 + ^||u - vf. 

By using m and the previous estimate, we derive the inequality 

\\\u-v\\\^ + f^\\Vy{u) - Vy{v)f - ^||u - vf 

F 

= \\\yiu) - y{v)f + -f^\\yy{u) - v?/(u)f + Q - Ik - vf 

< J{y{v),v) -J{y{u),u) + (||r- uV7?|| + Cf||/+ + divT||)^ . 

Finally, choosing <5 = 1/(2A) yields the estimate (IMll . □ 


Theorem 6. For any G V, let V(^ G Uad be given by the projection formula S21\) . Then, we 
obtain the following error majorant: 

(35) llk-U(;|||? < Mf{a,/3;r)X,T,p,Vi;) 

for arbitrary rj G V, t, p G H (div, O) and a, /3 > 0, where 

y 

(l + a)(l+/3)C|, 


A 4 f(a,/ 3 ;r?,C,'r,p,u<;) = ^\\r] - ydf + 


2a/3u^ 


11/+ + divrk - / {f + + divr) (dx 

Ja 


- [ (x - UV77) • VCda: + - (Cf Ik - yd+ divp|| + ||p-uVCII) 

Jn P 

X (^^Ik- yd + divp|| + ^||p- uVCII ACfIIZ + Uc; +divT|| + ||t- fV77||^ 


3A 


+ (Ik - + CfII/ + VC + div' 
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Proof. Applying (l34ll as well as (fTOl) and (l29l) finally leads to the estimate (l35l) . □ 

Proposition 2. The majorant At® (a, /3; p, defined in fi35\) attains the exact lower bound 

on the exact solution of the optimal control problem or, equivalently, of the optimality 

system i.e.. 


inf 

«C=P[u„.ui]{“<i-iC},a,;S >0 


Mf{a,P;r],C,T,p,V() = 0. 


The infimum is attained for V(^ = u, rj = y(u), f = p(u), t = vVyiu) and p = vVpiu). 

Proof. We have that 

Mfia,p-,y{u),p{u),vVy{u),vVp{u),u) = ^\\yiu) - yaf, 
which is zero if we let a go to zero. □ 


5. The Unconstrained Case 


In the unconstrained case, we have that C/ad = T^(U), i.e., vq satisfies the optimality condition 
(36) C + ~ =0 in 0, 

or, equivalently, vq = Ud — yC in The majorant dill) and minorant d27D simplify to 


(37) 


J®{a,(3-,r]X,T) :=i^||77-2/d||2 + - j/Viyf 


(1 + a)(l + ^)C|. II ^ 1 A , J- _l |2 , 1 iiaii 2 

||/ + Ud-^C + divr|| +^IICII, 


where J"® (a,/?; ??, C, t) = J"® (a,/?; ??, r, u^), and 


J^{V:C,r,p) = ^Wv-VdV + + _^ (j + Ud- ^C + divr^ fdx 

(38) +[ {T-iyVri)-VCdx--{CF\\v-yd + divp\\ + \\p-iyVC\\) 

Jn iL 

X - Vd + ddYp\\ + ^||p- i/VCII + CF||/ + Md- ^C + divx|| + ||x- , 

respectively. We obtain the following estimate using minorant (1551) : 

(39) J{y{u),u) > J^{r],C,T,p) Vt?, C £ U Vr, p e iJ(div, U). 


The minorant J7® is sharp, i.e., 

(40) sup J‘^{r],C,T,p) = J{y{u),u), 

7] ,P^ H {div ,Q) 

since the supremum of is attained for the optimal state ij = y and adjoint state C = p, and their 
corresponding exact fluxes r = and p = v'Vp. The optimal control is given by it = Ud — jp. 
In the unconstrained case, the proof of Theorem |3] provides not only an inequality but even an 
equation. This result is presented in the following theorem: 

Theorem 7. For any control function v £ L'^(fl), we have that 

(41) |||u- u||p =J{y{v),v) - J{y{u),u). 


Proof. The equation (HTl) follows by repeating the proof of Theorem [3] with applying equation ([HI) 
instead of inequality (|S]). □ 
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The norm ||| • ||| defined in (l30ll can be represented in terms of the state and the adjoint state 
(instead of the control) by using the optimality condition ([5]) for u and v as follows: 


(42) 


u-v\\\^ = ^\\y{u)-y{v)f+ ^\\u-vf 

= \\\y(.u) - y{v)f + ^\\ud - jpiu) - {ud - ^p(T))f 
= \Mu)-y{v)f + ^Mu)-p{v)\\\ 


Next, we present similar results to those of Theorem |4] and Proposition [T] for the unconstrained 
case. 


Theorem 8. Let Uad = L‘^{Ll). For any rj,(^ GV, we obtain the following error majorant: 


(43) 


, - UfllP < M®{a,l3;r)X,T,p) := J®{a,P;r],C,T) - J^{vX,T,p) 


with |||w — = \\\u — Ud jCWr 

M®{a,P;r],C,T,p) = ^\\v - ydf + - I'Vrjf 


(l + a)(l+/3)C|, 
2a/3u2 


11/ 


+ Ud-jC + divT\\'^ - J (^f + Ud- ^C + divr^ (dx 


- [ {t - uVry) -yCdxF- {Cp \\r] - yd + div p\\ + \\p- uVC||) 

Jq E 


n 
Cl 


Wv-yd + divpW + ^||p- uVCII +CF\\f + Ud - ic + divT|| + ||x- uVr^H ) , 


Au" ' Au 

where T,p G i^(div, fl) are arbitrary and a,/3 > 0. 


Proof. The estimate follows by applying (ED and then using the estimates and {|39j) . fH 


Proposition 3. Let Uad = Lf{Ll). The majorant {a, fd^rj^ ^ p) defined in attains the 
exact lower bound on the exact solution of the optimal control problem EP-iO), or, equivalently, of 
the optimality system and m, 


inf 

V,T,pG^/^(div,Q), 

a,^>0 


7W®(a,/3;77,C,'r, 


p)=0. 


In other words, the infimum is attained for the optimal state rj = y and adjoint state f = p with 
the optimal control u = Ud — jp, and for the exact fluxes t = uVy and p = uVp. 


Proof. We have that 


M®{a,l3-,y{u),p{u),v\Jy{u),v\Jp{u)) =^\\y{u) - ydW^, 
which is zero if we let a go to zero. 


□ 


Finally, we repeat the results of Theorems [S] and Proposition [5] for the unconstrained case. 


Theorem 9. Let Uad = Lf{Ll). For any gV , we obtain the following error majorant: 
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with |||u — = \\\u — Ud + jCllli 

M^{a,l3-,ri,C„T,p) = - Vdf + t H'^ ~ 


2q;i/2 


(1 + a;)(l + (i)Cp 
2aj3i^ 


\\f + Ud - ic + divxf - J (^f + Ud- ^C + divT^ (dx 


- [ (t - ^V77) • VCda: + - (Cf Ih-?/d+ divp|| + ||p-i/VCII) 
Jn IL 


Cl 


Cl 


X ( -^h - yd + divp|| + ^||p- i^VCII +CF\\f + Ud - -C + divT|| + ||t - FVy|| 


+ 


3A 


\t - vVriW + CfII/ + Md - + divrl 


for arbitrary t, p G H (div, O) and a, P > 0. 


Proof. Applying ((Ml) in the unconstrained case as well as (flUl) and (IMl) finally leads to the estimate 

dm). □ 


Proposition 4. Let Uad = {Vi). The majorant Mf {a, j3\rj^C,^T^ p) defined in attains the 
exact lower bound on the exact solution of the optimal control problem (QP-®, or, equivalently, of 
the optimality system and m, *-e., 


inf 

7] ,(^^V,T, pG H {div .Q), 
a,l3>0 


Mf{a,P;r]X,T, 


p)=0. 


The infimum is attained for the optimal state tj = y and adjoint state C = p with the optimal 
control u = Ud — jP, and for the exact fluxes T = vVy and p = FVp. 


Proof. We have that 

Mf {a, l3;y{u),p{u),vVy{u),vVp{u)) = ^ ||y(M) - y^f, 
which is zero if we let a go to zero. □ 


6. The Finite Element Discretization and the Preconditioned MINRES Solver 


In this section, we present the finite element discretization of the optimality system in order to 
derive approximations of the state, control and adjoint state. These approximations can be used 
for the computation of the majorants and minorants. We start with the unconstrained case. 


The unconstrained case. Since the control can be eliminated from the optimality system by 
using m, we only have to solve the system (ED-© for the state y GV and the adjoint state p G V. 
For that, we approximate these unknown functions by hnite element functions yh, Ph G Vh G V, 
where the finite element space Vh is defined as 

(45) Vd = span{ 


with the standard nodal basis {ipi{x) = ipihix) : i = 1, 2,..., Uh} and h denotes the discretization 
parameter (mesh size) such that n = Uh = dimW = 0{h~^). Using continuous, piecewise linear 
finite elements on triangles on a regular triangulation to construct the hnite element subspace Vh 
and its basis (see, e.g., H]), yields the following linear system: 


(46) 


Mh -Kh 
-Kh -X-^Mh 


Ph 


.h - 


Pd 


-if+ Ud) 


which has to be solved with respect to the nodal parameter vectors y^ = (yi)i=i,,,,,n G R" and 


p^ = (pi)i=i,...,n G M" of the hnite element approximations 

n n 

yhix) = '^yi(pflx) and Phjx) = y^pjipflx) 


2=1 


2=1 


(47) 
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to the unknown functions y{x) and p{x). The matrices Mh and Kh correspond to the mass and 
stiffness matrices, respectively. Their entries are computed by the following formulas: 






V Vtpi • 'S/ipj dx, 


i,j = 


The right hand sides are given by 


y^= yd dx 


Uw = 


'-JQ. 


Ud y>j dx 




and / = 




fipj dx 


J j=l,...,n 


The discretized optimality system (l46ll is a saddle point problem and can be solved by a precon¬ 
ditioned MINRES (minimal residual) method, see [55] . Hence, it is crucial to construct precondi¬ 
tioners, which yield robust and fast convergence for the preconditioned MINRES method. In EZ], 
the following preconditioner was constructed: 


(48) 


V = 


Mh + VXKh 

0 


iMh + ^^Kh 


leading to a robust convergence with respect to h and A (as well as v for our problem). 


The constrained case. In the case of having inequality constraints imposed on the control, we 
can reformulate the variational inequality (I5|), which is equivalent to the projection formula 

uix^ ^[Ua{x),Ub{x)]\'^d{x') —p(3l)}, 

by introducing an additional parameter p e see, e.g., [T51155] . The idea is to apply a 

primal-dual active set strategy in order to linearize the optimality system, see, e.g., [T5]. The 
parameter p is defined as /r = —jp + Ud — u for the exact solution u € C/ad- Moreover, if a function 
V € C/ad satisfies the relations 




Ua{x) 

if v{x) -1- pl{x) < Ua{x), 


v{x) = < 

Ud{x) - Xpix) 

if v{x) -1- p{x) G [Ua{x), Ubix)], 



Ub{x) 

if vix) + p{x) > Ubix), 

then it is the optimal solution, since this means that v satisfies the projection formula. In the 
following, we want to state the main steps of the primal-dual active set method for the optimal 
control problem ©-([SI). Given an iterate (yh-i,Uh-i,Pk-i, Pk-i), the active and inactive sets are 


determined as follows: 

Ah = {x € fl : Uk(x) -I- fXk(x) < Ua(x)}, Ai = {x Gil : Uk{x) + Pkix) > Ub{x)}, 

ik = n\iAtuAi). 


If A% = A%_i and Ah = Ah-i, then we have attained the optimal solution. Otherwise, the next 
iterate is the solution of the system 

-div{i3{x)\7yk{x)) = f{x) + Uk{x), x G il, yk{x)=0,xGT, 

-div{i3{x)\7pk{x)) = yh{x) - yd{x), x G il, ph{x) = 0, a: £ E, 

Uh{x) - Udix) + X~^XikPk{x) = XAiUaix) + x^6Ub(a;), xGil, 

where x denotes the characteristic function. We eliminate again the control Uh from the system 
and obtain the following reduced optimality system written in its variational formulation: 


/ yhwdx— / ph ■ Vw dx = / ydwdx, 

/ vS/yh-S/w dx / phwdx= {f + Ud)w dx + / Uawdx+ Ubwdx, 
Jn Jik Jn JAf JAl 
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for all w,q € V. Now, discretizing the reduced optimality system by the finite element method 
with the space Vh <ZV defined as in (H51) leads to the linear system 


(49) 


( Mh -Kh \ 

\ -Kh -X-^Mh,i, ) 


where the matrix M^^Xk with the entries 


Ph,k 




-(/ + Ud) -Ua - Wfe 


and the vectors 


= 


’Ik 




i dx, 


= 1,... ,n, 


Ua = 

/ Ua <fj dx 

and Uh = 

/ Ub <Pj dx 




Mai J 


have to be computed in every fc-th iteration step of the level set method. In order to solve the 
linear system by the preconditioned MINRES method, we need again a preconditioner, which 
provides fast convergence and is as robust as possible. Here, we refer to the preconditioner 


(50) 


-Px, = 


Mh + '/XKh 
0 


jMh,Xk 




which was presented in m- For the preconditioner (1501) , it can be proven robustness with respect 
to the inactive set Ik and the mesh size h (as well as with respect to the diffusion parameter v in 
our case), but not with respect to the cost parameter A (the condition number scales like l/\/X). 


7 . Numerical Results 

In this section, we present and discuss first numerical results for the unconstrained case. The 
computational domain O = (0,1) x (0,1) is uniformly decomposed into triangles, and standard con¬ 
tinuous, piecewise linear finite elements are used for the discretization. In this case, the Friedrichs 
constant is Cp = l/('\/27i'). 

The construction of rj, t and p is an important issue in order to obtain sharp guaranteed 
bounds from the majorants and minorants. Since rj and C are constructed by continuous, piecewise 
linear approximations, their gradients are only piecewise constant. Then, V?7, VC € L^(H), but 
V?7, VC ^ II(div,fi). Hence, a flux reconstruction is needed in order to obtain suitable fluxes 
T, p G H{div, H). A good reconstruction of the flux is an important and nontrivial topic. We can 
regularize the fluxes by a post-processing operator which maps the L^-functions into H(div, H), 
see [3S]. There are various techniques for realizing these post-processing steps such as, e.g., local 
post-processing by an elementwise averaging procedure or by using Raviart-Thomas elements, see 
[3S1I25] and references, therein. In our numerical experiments, we use Raviart-Thomas elements 
of the lowest order, see, e.g., [29l[5l|36]. We define the normal fluxes on interior edges Emn by 

(t ■ = (Amn(Vr7)|T„ + (1 - Xmn){yil)\Tj ' 

(p ■ nE,nn)\E^„ = (Amn(VC)|T„ + (1 “ Amn)(VC)|T„) ' i 

with Amn = 1/2 due to uniform discretization. Here, (V?7 )|t„, (VC)|rmj (V77 )|t„ and (VC)|r„ 
are constant vectors on two arbitrary, neighboring elements Tm and r„. On bonndary edges, the 
only one existing flux is used. Hence, three normal fluxes are defined on the three sides of each 
element. Inside, we reconstruct the fluxes r and p by the standard lowest-order Raviart-Thomas 
(RT°-) extension of normal fluxes with 

RT°(r^) :={Te (T2 (t)) 2 :VTe7/ 3a,6,cGM Va; G T, 

t{x) = (a, bX + cx and [t]e • ue = 0 V interior edges E}, 

where [t]^; denotes the jump of r across the edge E shared by two neighboring elements on a 
triangulation Th- Altogether, the RT°-extension yields averaged fluxes from H(div, H), i.e., 

T = GRT(Vr 7 ), p = Grt(VC), Grt : T^(H) —?> H(div, H). 
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Results for functional a posteriori error estimates of parabolic time-periodic boundary value prob¬ 
lems using the same discretization and flux reconstruction techniques can be found in mi- 

In order to solve the saddle point systems (H51) . we use the AMLI preconditioner proposed by 
Kraus in [TB] for an inexact realization of the block-diagonal preconditioner (H51) in the MINRES 
method. The numerical results where computed on grids of different mesh sizes (from 8x8 
to 256 X 256). The preconditioned MINRES iteration was stopped after 8 iteration steps in all 
computations using the AMLI preconditioner with 4 inner iterations. The presented CPU times in 
seconds include the computational times for computing the majorants and minorants, which 
are very small in comparison to the computational times of the solver. All computations were 
performed on a laptop with Intel(R) Core(TM) i5-4308U CPU @ 2.80GHz. 

In the numerical experiment, we consider the following given data: the desired state 

yd{x) = (1 + 0.0471^) sin(a;i7r) sin(a;27r), 


the desired control Ud{x) = 0 and the source term f{x) = 0. We choose the parameter v = \ and 
the cost parameter A = 0.01. The exact state and control are known and given by 

y{x) = sin(a;i7r) sin(x27r) and u{x) = 2?!^ sin(aii7r) sin(a:27r), 


respectively. Hence, the exact value of the cost functional is given by J{y{u)^u) = 2.385. Table[T] 
presents the CPU times in seconds U®®, the majorants ([H7| and minorants ([BB]! as well as the 
efficiency indices 


J{y{u),u) 

J®/Q ^ 

®ff je{yX,T,p) 


I® - 

-'eff — 


J{y{u):U) 


tMi _ 
^eff — 


I Mf{aX;r]X,x) 


|M-'ycllli 


obtained on grids of different mesh sizes. The parameters a and /3 are chosen by a proper mini¬ 
mization of the majorant J® with respect to a and (3, see m- In Table [BJ the individual parts 
of the majorants and minorants are presented. They are denoted as follows: 


= lldivr-icil, T^i = \\v-yd + <iivp\\, = ||r - VryH, 7^^ = ||p-VC||, 

T^ = \\v-yd\\, 7^^ = ||CII, TZ 4 = J^(^divT - Cdx, = J^{t - vVvi)-VQdx. 


grid 

U®® 


J® 

T® 

-'eff 

T® 

-'eff 

7-®/e 

-'eff 

tMi 

-'eff 

8 

X 

8 

0.006 

2.248 

2.438 

1.022 

1.061 

1.085 

2.369 

16 

X 

16 

0.012 

2.351 

2.431 

1.019 

1.015 

1.034 

1.953 

32 

X 

32 

0.045 

2.376 

2.412 

1.011 

1.004 

1.015 

1.751 

64 

X 

64 

0.179 

2.383 

2.399 

1.006 

1.001 

1.007 

1.656 

128 

X 

128 

0.784 

2.384 

2.392 

1.003 

1.000 

1.003 

1.610 

256 

X 

256 

3.246 

2.385 

2.389 

1.002 

1.000 

1.002 

1.587 


Table 1. Efficiency of the minorants and majorants. 


In both tables, we observe the efficiency of the AMLI preconditioned MINRES method presented 
in [HHD]. The computational times increase with a factor of four. Moreover, one can see that 
Tidl and 'R!{ reduce as a factor of two, and iVl and as a factor of four showing the efficiency of 
the applied flux reconstruction. However, one could consider other flux reconstruction techniques 
in terms of improving the efficiency indices such as using higher order Raviart-Thomas elements, 
see Altogether we can observe that the majorants and minorants provide good estimates for 
the value of the cost functional and can be used in order to compute guaranteed estimates for the 
discretization error in the combined norms as it is discussed and proved in this paper. 
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grid 


n\ 


t4 

^3 

^3 

77-4 

77-5 

8x8 

1.146 

0.235 

0.053 

0.011 

1.928 

0.094 

-0.013 

0.001 

16 X 16 

0.606 

0.121 

0.016 

0.003 

1.943 

0.098 

-0.004 

0.000 

32 X 32 

0.307 

0.061 

0.004 

0.001 

1.947 

0.098 

-0.001 

0.000 

64 X 64 

0.154 

0.030 

0.001 

0.000 

1.948 

0.099 

-0.000 

0.000 

128 X 128 

0.077 

0.015 

0.000 

0.000 

1.948 

0.099 

-0.000 

0.000 

256 X 256 

0.039 

0.008 

0.000 

0.000 

1.948 

0.099 

-0.000 

0.000 


Table 2. The indiviudal parts of the minorants and majorants. 


8. Conclusions 

This work was devoted to the derivation of guaranteed and fully computable lower bounds 
(minorants) for cost functionals of distributed elliptic optimal control problems in order to close 
the gap of the already existing results on upper bounds for these cost functionals, see [?]• An 
important result of this work was to prove that the discretization error in the state and the 
control can be computed by the difference between majorant and minorant of the cost functional. 
Altogether we derive a fully computable upper bound for the discretization error in the state and 
the control, which can be, in principle, used as object of direct minimization. However, the aim 
of this work was not to present an algorithm for this minimization, but to derive first results 
on minorants for cost functionals of distributed elliptic optimal control problems with control 
constraints as well as present first numerical tests for the theoretical results derived. 
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